skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Search for: All records

Creators/Authors contains: "Bellovary, Jillian"

Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher. Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?

Some links on this page may take you to non-federal websites. Their policies may differ from this site.

  1. Abstract I hypothesize a physical explanation for the “little red dots” (LRDs) discovered by the James Webb Space Telescope (JWST). The first star formation in the Universe occurs in dense clusters, some of which may undergo runaway collapse and form an intermediate mass black hole. This process would appear as a very dense stellar system, with recurring tidal disruption events (TDEs) as stellar material is accreted by the black hole. Such a system would be compact, UV-emitting, and exhibit broad Hαemission. If runaway collapse is the primary mechanism for forming massive black hole seeds, this process could be fairly common and explain the large volume densities of LRDs. In order to match the predicted number density of runaway collapse clusters, the tidal disruption rate must be on the order of 10−4per year. A top-heavy stellar initial mass function may be required to match observations without exceeding the predicted ΛCDM mass function. The TDE LRD hypothesis can be verified with follow-up JWST observations looking for TDE-like variability. 
    more » « less
    Free, publicly-accessible full text available May 6, 2026
  2. Abstract A consequence of a nonzero occupation fraction of massive black holes (MBHs) in dwarf galaxies is that these MBHs can become residents of larger galaxy halos via hierarchical merging and tidal stripping. Depending on the parameters of their orbits and original hosts, some of these MBHs will merge with the central supermassive black hole in the larger galaxy. We examine four cosmological zoom-in simulations of Milky Way-like galaxies to study the demographics of the black hole mergers that originate from dwarf galaxies. Approximately half of these mergers have mass ratios less than 0.04, which we categorize as intermediate mass ratio inspirals, or IMRIs. Inspiral durations range from 0.5–8 Gyr, depending on the compactness of the dwarf galaxy. Approximately half of the inspirals may become more circular with time, while the eccentricity of the remainder does not evolve. Overall, IMRIs in Milky Way-like galaxies are a significant class of black hole mergers that can be detected by LISA, and must be prioritized for waveform modeling. 
    more » « less
    Free, publicly-accessible full text available June 19, 2026
  3. Abstract We examine the quenching characteristics of 328 isolated dwarf galaxies 10 8 < M star / M < 10 10 within theRomulus25cosmological hydrodynamic simulation. Using mock-observation methods, we identify isolated dwarf galaxies with quenched star formation and make direct comparisons to the quenched fraction in the NASA Sloan Atlas (NSA). Similar to other cosmological simulations, we find a population of quenched, isolated dwarf galaxies belowMstar< 109Mnot detected within the NSA. We find that the presence of massive black holes (MBHs) inRomulus25is largely responsible for the quenched, isolated dwarfs, while isolated dwarfs without an MBH are consistent with quiescent fractions observed in the field. Quenching occurs betweenz= 0.5–1, during which the available supply of star-forming gas is heated or evacuated by MBH feedback. Mergers or interactions seem to play little to no role in triggering the MBH feedback. At low stellar masses,Mstar≲ 109.3M, quenching proceeds across several Gyr as the MBH slowly heats up gas in the central regions. At higher stellar masses,Mstar≳ 109.3M, quenching occurs rapidly within 1 Gyr, with the MBH evacuating gas from the central few kpc of the galaxy and driving it to the outskirts of the halo. Our results indicate the possibility of substantial star formation suppression via MBH feedback within dwarf galaxies in the field. On the other hand, the apparent overquenching of dwarf galaxies due to MBH suggests that higher-resolution and/or better modeling is required for MBHs in dwarfs, and quenched fractions offer the opportunity to constrain current models. 
    more » « less
  4. Abstract We explore the characteristics of actively accreting massive black holes (MBHs) within dwarf galaxies in the Romulus25cosmological hydrodynamic simulation. We examine the MBH occupation fraction, X-ray active fractions, and active galactic nucleus (AGN) scaling relations within dwarf galaxies of stellar mass 108M<Mstar< 1010Mout to redshiftz= 2. In the local universe, the MBH occupation fraction is consistent with observed constraints, dropping below unity atMstar< 3 × 1010M,M200< 3 × 1011M. Local dwarf AGN in Romulus25follow observed scaling relations between AGN X-ray luminosity, stellar mass, and star formation rate, though they exhibit slightly higher active fractions and number densities than comparable X-ray observations. Sincez= 2, the MBH occupation fraction has decreased, the population of dwarf AGN has become overall less luminous, and as a result the overall number density of dwarf AGN has diminished. We predict the existence of a large population of MBHs in the local universe with low X-ray luminosities and high contamination from X-ray binaries and the hot interstellar medium that are undetectable by current X-ray surveys. These hidden MBHs make up 76% of all MBHs in local dwarf galaxies and include many MBHs that are undermassive relative to their host galaxy’s stellar mass. Their detection relies on not only greater instrument sensitivity but also better modeling of X-ray contaminants or multiwavelength surveys. Our results indicate that dwarf AGN were substantially more active in the past, despite having low luminosity today, and that future deep X-ray surveys may uncover many hidden MBHs in dwarf galaxies out to at leastz= 2. 
    more » « less
  5. ABSTRACT Stars and stellar remnants orbiting a supermassive black hole (SMBH) can interact with an active galactic nucleus (AGN) disc. Over time, prograde orbiters (inclination i < 90°) decrease inclination, as well as semimajor axis (a) and eccentricity (e) until orbital alignment with the gas disc (‘disc capture’). Captured stellar-origin black holes (sBH) add to the embedded AGN population that drives sBH–sBH mergers detectable in gravitational waves using LIGO–Virgo–KAGRA or sBH–SMBH mergers detectable with Laser Interferometer Space Antenna. Captured stars can be tidally disrupted by sBH or the SMBH or rapidly grow into massive ‘immortal’ stars. Here, we investigate the behaviour of polar and retrograde orbiters (i ≥ 90°) interacting with the disc. We show that retrograde stars are captured faster than prograde stars, flip to prograde orientation (i < 90°) during capture, and decrease a dramatically towards the SMBH. For sBH, we find a critical angle iret ∼ 113°, below which retrograde sBH decay towards embedded prograde orbits (i → 0°), while for io > iret sBH decay towards embedded retrograde orbits (i → 180°). sBH near polar orbits (i ∼ 90°) and stars on nearly embedded retrograde orbits (i ∼ 180°) show the greatest decreases in a. Whether a star is captured by the disc within an AGN lifetime depends primarily on disc density, and secondarily on stellar type and initial a. For sBH, disc capture time is longest for polar orbits, low-mass sBH, and lower density discs. Larger mass sBH should typically spend more time in AGN discs, with implications for the spin distribution of embedded sBH. 
    more » « less
  6. ABSTRACT We present a phenomenological forward Monte Carlo model for forecasting the population of active galactic nuclei (AGNs) in dwarf galaxies observable via their optical variability. Our model accounts for expected changes in the spectral energy distribution of AGNs in the intermediate-mass black hole (IMBH) mass range and uses observational constraints on optical variability as a function of black hole (BH) mass to generate mock light curves. Adopting several different models for the BH occupation function, including one for off-nuclear IMBHs, we quantify differences in the predicted local AGN mass and luminosity functions in dwarf galaxies. As a result, we are able to model the fraction of variable AGNs as a function of important galaxy host properties, such as host galaxy stellar mass, in the presence of selection effects. We find that our adopted occupation fractions for the ‘heavy’ and ‘light’ initial BH seeding scenarios can be distinguished with variability at the 2–3σ level for galaxy host stellar masses below ∼108M⊙ with data from the upcoming Vera C. Rubin Observatory. We also demonstrate the prevalence of a selection bias whereby recovered IMBH masses fall, on average, above the predicted value from the local host galaxy–BH mass scaling relation with the strength of this bias dependent on the survey sensitivity. Our methodology can be used more broadly to calibrate AGN demographic studies in synoptic surveys. Finally, we show that a targeted ∼ hourly cadence program over a few nights with the Rubin Observatory can provide strong constraints on IMBH masses given their expected rapid variability time-scales. 
    more » « less
  7. Abstract Current observational facilities have yet to conclusively detect 10 3 –10 4 M ⊙ intermediate-mass black holes (IMBHs) that fill in the evolutionary gap between seed black holes in the early universe and z ∼ 0 supermassive black holes. Dwarf galaxies present an opportunity to reveal active IMBHs amidst persistent star formation. We introduce photoionization simulations tailored to address key physical uncertainties: coincident versus noncoincident mixing of IMBH and starlight excitation, open versus closed geometries of surrounding gas clouds, and different shapes of the spectral energy distribution of active galactic nuclei (AGN). We examine possible AGN emission line diagnostics in the optical and mid-IR, and find that the diagnostics are often degenerate with respect to the investigated physical uncertainties. In spite of these setbacks, and in contrast to recent work, we are able to show that [O iii ]/H β typically remains bright for dwarf AGN powered by IMBHs down to 10 3 M ⊙ . Dwarf AGN are predicted to have inconsistent star-forming and Seyfert/LINER classifications using the most common optical diagnostics. In the mid-IR, [O iv ] 25.9 μ m and [Ar ii ] 6.98 μ m are less sensitive to physical uncertainties than are optical diagnostics. Based on these emission lines, we provide several diagrams of mid-IR emission line diagnostic diagrams with demarcations for separating starbursts and AGN with varying levels of activity. The diagrams are valid over a wide range of ionization parameters and metallicities out to z ∼ 0.1, so will prove useful for future JWST observations of local dwarf AGN in the search for IMBHs. We make our photoionization simulation suite freely available. 
    more » « less
  8. Abstract Existing star-forming vs. active galactic nucleus (AGN) classification schemes using optical emission-line diagnostics mostly fail for low-metallicity and/or highly star-forming galaxies, missing AGN in typicalz∼ 0 dwarfs. To recover AGN in dwarfs with strong emission lines (SELs), we present a classification scheme optimizing the use of existing optical diagnostics. We use Sloan Digital Sky Survey emission-line catalogs overlapping the volume- and mass-limited REsolved Spectroscopy Of a Local VolumE (RESOLVE) and Environmental COntex (ECO) surveys to determine the AGN percentage in SEL dwarfs. Our photoionization grids show that the [Oiii]/Hβversus [Sii]/Hαdiagram (Siiplot) and [Oiii]/Hβversus [Oi]/Hαdiagram (Oiplot) are less metallicity sensitive and more successful in identifying dwarf AGN than the popular [Oiii]/Hβversus [Nii]/Hαdiagnostic (Niiplot or “BPT diagram”). We identify a new category of “star-forming AGN” (SF-AGN) classified as star-forming by the Niiplot but as AGN by the Siiand/or Oiplots. Including SF-AGN, we find thez∼ 0 AGN percentage in dwarfs with SELs to be ∼3%–16%, far exceeding most previous optical estimates (∼1%). The large range in our dwarf AGN percentage reflects differences in spectral fitting methodologies between catalogs. The highly complete nature of RESOLVE and ECO allows us to normalize strong emission-line galaxy statistics to the full galaxy population, reducing the dwarf AGN percentage to ∼0.6%–3.0%. The newly identified SF-AGN are mostly gas-rich dwarfs with halo mass <1011.5M, where highly efficient cosmic gas accretion is expected. Almost all SF-AGN also have low metallicities (Z≲ 0.4Z), demonstrating the advantage of our method. 
    more » « less
  9. null (Ed.)
    ABSTRACT Active galactic nuclei (AGN) are powered by the accretion of discs of gas on to supermassive black holes (SMBHs). Stars and stellar remnants orbiting the SMBH in the nuclear star cluster (NSC) will interact with the AGN disc. Orbiters plunging through the disc experience a drag force and, through repeated passage, can have their orbits captured by the disc. A population of embedded objects in AGN discs may be a significant source of binary black hole mergers, supernovae, tidal disruption events, and embedded gamma-ray bursts. For two representative AGN disc models, we use geometric drag and Bondi–Hoyle–Littleton drag to determine the time to capture for stars and stellar remnants. We assume a range of initial inclination angles and semimajor axes for circular Keplerian prograde orbiters. Capture time strongly depends on the density and aspect ratio of the chosen disc model, the relative velocity of the stellar object with respect to the disc, and the AGN lifetime. We expect that for an AGN disc density $$\rho \gtrsim 10^{-11}{\rm g\, cm^{-3}}$$ and disc lifetime ≥1 Myr, there is a significant population of embedded stellar objects, which can fuel mergers detectable in gravitational waves with LIGO-Virgo and LISA. 
    more » « less
  10. null (Ed.)
    ABSTRACT Massive black holes often exist within dwarf galaxies, and both simulations and observations have shown that a substantial fraction of these may be off-centre with respect to their hosts. We trace the evolution of off-centre massive black holes (MBHs) in dwarf galaxies using cosmological hydrodynamical simulations, and show that the reason for off-centre locations is mainly due to galaxy–galaxy mergers. We calculate dynamical time-scales and show that off-centre MBHs are unlikely to sink to their galaxys’ centres within a Hubble time, due to the shape of the hosts’ potential wells and low stellar densities. These wandering MBHs are unlikely to be detected electromagnetically, nor is there a measurable dynamical effect on the galaxy’s stellar population. We conclude that off-centre MBHs may be common in dwarfs, especially if the mass of the MBH is small or the stellar mass of the host galaxy is large. However, detecting them is extremely challenging, because their accretion luminosities are very low and they do not measurably alter the dynamics of their host galaxies. 
    more » « less